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The hole-doped ground state of a recently introduced extended "t-U-J" two-orbital Hubbard 
model for the Fe-based superconductors is studied via exact diagonalization methods on small clus- 
ters. Similarly as in the previously studied case of electron doping, A. Nicholson et at, Phys. Rev. 
Lett. 106 21702 (2011), upon hole doping it is observed that there are several competing pairing 
symmetries including Aig, Big, and B2g. However, contrary to the electron-doped case, the ground 
state of the hole-doped state has pseudocrystal momentum k = {n,n) in the unfolded Brillouin 
zone. In the two Fe-atom per unit cell representation, this indicates that the ground state involves 
anti-bonding, rather than bonding, combinations of the orbitals of the two Fe atoms in the unit-cell. 
The lowest state with k = (0, 0) has only a slightly higher energy. These results indicate that this 
simple two-orbital model may be useful to capture some subtle aspects of the hole-doped pnictides 
since calculations for the five-orbital model have unveiled a hole pocket centered at M (k = (7r,7r)) 
in the unfolded Brillouin zone. 

PACS numbers: 74.20.Rp, 71.10.Fd, 74.70.Xa, 75.10.Lp 



I. INTRODUCTION 

The detailed study of the recently discovered high crit- 
ical temperature superconductivity (HTCS) in the iron- 
based pnictides and chalcogenides^ continues providing 
important information to understand the still puzzling 
mechanism that drives this remarkable phenomenon. Ex- 
periments indicate that these Fe-based materials share 
many properties with the high cuprateS)^ such as mag- 
netically ordered parent compounds^ and superconduct- 
ing states stabilized by either electron or hole dopingii 
However, there are also remarkable differences, such as 
the fact that the parent compounds, at least for the 
pnictides, are (bad) metals rather than Mott insulators. 
Moreover, several of the iron d-orbitals are active at 
the Fermi surface (FS) as opposed to the case of the 
cuprates where just the copper dx2_y2 orbital plays the 
major role. While it has been clearly established that 
the superconducting pairing operator has d-wave sym- 
metry in the case of the hole-doped cuprates?^ the sym- 
metry of the pairing operator in the pnictides is still con- 
troversial: surface-sensitive angle-resolved photoemission 
(ARPES) studies^ indicate that full nearly-momentum- 
independcnt gaps open on all FS pockets, compatible 
with the s± statci^ However, several other experiments 
testing bulk properties provide results compatible with 
nodal superconductivity.— 

Reliable theoretical studies are difficult to implement 
for the complex multi-orbital models needed for these 
materials without making explicit assumptions about the 
ground state properties or about the mechanism and 
strength of the pairing interactions. Under the assump- 
tion of a magnetically driven superconducting instabil- 
ity, the Random Phase Approximation (RPA) is often 
applied in this context providing indications that several 



pairing channels, mostly Aig and -Big, are in competi- 
tion in this type of compounds 4 However, RPA relies on 
a particular subset of Feynman diagrams and it is a weak- 
coupling approach. On the opposite extreme, strong cou- 
pling studies have also been performed. Depending on 
the particular model used, some authors have found ev- 
idence of pairing with i?2g symmetr}*^ while others have 
found a variety of competing statesJ^ A complementary 
approach to these previous investigations is to perform 
an exact diagonalization (ED) of the model Hamiltoni- 
ans, allowing to solve the problem exactly for any value of 
the interaction. However, since the Hilbert space grows 
exponentially with the system size, this method can be 
implemented only in very small clusters and with a re- 
duced number of active Fe orbitals. This ED approach 
has been recently applied by the authors and collabora- 
tors to the study of an electron-doped two-orbital Hub- 
bard model and in that effort the presence of competing 
pairing states was observed as the strength of the cou- 
pling parameters was varied.—"— 

An important characteristic of the widely studied Hub- 
bard models for the pnictides/chalcogenides is that they 
are not particle-hole symmetric. On the experimental 
side, superconductivity has been found both upon elec- 
tron and hole doping, but it seems that hole-doped mate- 
rials belonging to the 122 family are more suitable for the 
use of surface-sensitive techniques, while electron-doped 
materials belonging to the 1111 family are more easily 
studied with bulk techniques i^ii Then, it is natural to 
wonder whether a potential source of the differences in 
the experimental results regarding the pairing symme- 
tries may arise from the nature of the dopants. For this 
reason it is important to study theoretically the prop- 
erties of multi-orbital Hubbard models both under elec- 
tron and hole doping. Previously, such a comparative 
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analysis has been performed employing the RPA method 
applied to a five-orbital Hubbard modelJ^ For the pa- 
rameter range studied in that case (weak coupling), a 
pairing state with Aig symmetry was observed in both 
cases. The state found has nodes on the electron pockets 
in the electron-doped case, but upon hole-doping an extra 
hole pocket at the M = (tt, tt) point in the Brillouin zone 
leads to the removal of the nodes and the development 
of a nodeless s± stateJ^ 

The goal of the present publication is to study the 
most favorable pairing channels of a two-orbital Hubbard 
model using small-cluster exact diagonalization tech- 
niques (namely, the Lanczos algorithm) for the case of 
hole doping, and to contrast the results against those 
found for the case of electron-doping of the same model 
that have been recently reported.—"— To reduce the se- 
vere constraints imposed by the small size of the clusters 
that can be diagonalized in present day computers, a sim- 
ple generalization of the Hubbard model for the pnictides, 
previously introduced for the study of the electron doped 
case^^ will be here applied. For this purpose, Heisenberg 
"J" terms will be added to the original Hubbard model to 
enhance spin order and pairing tendencies, but without 
projecting out doubly occupied sites and charge fluctua- 
tions. These terms help to establish tightly bound-states 
upon doping that can be studied with Lanczos methods 
on the small clusters currently accessible with state-of- 
the-art computers. Comparisons with RPA results for 
the five-orbital model and with experimental data will 
be performed. 

The organization of the paper is as follows: the model 
and the method are reviewed in Section [Hi the main re- 
sults for hole doping are presented in section IHIl while 
section HVl is devoted to the conclusions. 



II. MODEL AND METHOD 

The model studied here is based on the well-known 
and widely used two-orbital Hubbard mode l^^i^^'^^ that 
employs the dxz (x) and dyz (y) Fe orbitals. These or- 
bitals provide the largest contribution to the pnictides' 
band structure at the FS.— The reduction in the actual 
number of active orbitals in the pnictides is necessary in 
order to perform the present Lanczos studies. Calcula- 
tions with more orbitals for the same cluster studied here 
are simply not possible at present. 

The parameters of the electronic hopping terms of the 
model were previously chosen to provide a close agree- 
ment with the band structure calculations obtained with 
density-functional theoryJ^ In addition to the hopping 
terms, the model also includes the on-site Coulomb in- 
teraction consisting of intra- and inter-orbital Coulomb 
repulsions with couplings U and U' , the Hund's rule cou- 
pling Jh, and the pair-hopping term with strength J'. 
While U can in principle depend on the particular or- 
bital due to different screening effects for each orbital, 
this is not the case for the d^z and dyz orbitals that 



form a degenerate Cg doublet, implying that the rela- 
tions U'=U — 2J^, and J'= Jh are fulfilled for symme- 
try reasons.— This Hamiltonian, with the kinetic-energy 
hopping and onsite-interaction terms above mentioned, 
has been studied in detail previously! ^^•^'^i-'^^'^^ Moreover, 
in a recent investigation of the electron-doped case^^i the 
model was supplemented by a Heisenberg interactions to 
amplify the strength of the magnetic state and, conse- 
quently, the pairing strength as well. 

Naively, it may seem that selecting a stronger on-site 
Hubbard interaction would stabilize a stronger antiferro- 
magnetic state. However, this procedure also induces an 
insulator, and actually the strength of the effective cou- 
pling between the Fe-spins decreases as l/U with increas- 
ing U. To avoid this problem, in early studies of the one- 
band t-U-J modcii^ and in a recent study of the electron- 
doped two-orbital model,— Heisenberg terms have been 
added and shown to enhance pairing tendencies. Since 
our aim is to investigate the symmetries of the Cooper 
pairs, the additional magnetic interactions must have the 
same symmetries as the original Hamiltonian. To make 
sure that the symmetries are properly handled, the ad- 
ditional Heisenberg interaction is given by the operato- 
rial form that corresponds to the superexchange terms 
derived from the strong-coupling (large-?/) limit. In the 
case of the one-band Hubbard model, this is a Heisenberg 
term with spin S ~ 1/2. In the case of a multi-orbital 
model away from half-filling, the corresponding superex- 
change contains an orbital degree of freedom in addition 
to the spin and it is of a Kugcl-Khomskii typcj^i^i 

In the present case of a half-filled two-orbital model, 
the low-energy Hilbert space for the strong-coupling 
limit, with both U and Jh large, is given by doubly oc- 
cupied sites with singly occupied orbitals. Due to the 
Hund's coupling, the two electrons per site form a triplet 
state, with an energy Eq = U' — Jn = U — 3Jh, com- 
pared to El ^ U' + Ju = U - J' = Eq + 2Jh and 
E2 = U + Ju ^ Eq + 4 Jh for inter- and intra-orbital sin- 
glet states. The low-energy Hilbert space is, thus, given 
by a spin S* = 1 at each site. The interaction between 
these spins can be obtained by second-order perturba- 
tion theory in an analogous manner as the well-known 
derivation of the Heisenberg model from the one-orbital 
Hubbard model. The calculation for the two orbitals is 
the most easily carried out when the hopping term pre- 
serves orbital flavor, because the flrst hopping process, 
which creates a virtual excitation with energy U + Jn, 
then has to involve the same orbital as the second, which 
goes back to the low-energy Hilbert subspace. By this 
procedure it can be shown that the result is the isotropic 
Heisenberg interaction for S = I with a coupling 



where ta and tt are the hopping parameters correspond- 
ing to the two orbitals. With the notation tjj/f, = ti/2r^ 
the nearest-neighbor (NN) coupling Jnn can be derived. 
For a next-ncarest-neighbor (NNN) coupling, which is 
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natural since in the original Hubbard model the hoppings 
involve both NN and NNN Fe atoms, it is convenient to 
transform to a rotated orbital basis {\xz) ± \yz))/^/2^ 
where the hoppings are diagonal in orbital space and 
given by ± <4 , leading to 



>^NNN — -T^- — = ^ 



3 [/ + Jh t{ + t 



2 ^NN 



(2) 



Similarly as in the previous investigation of Ref. [llj 
and to avoid the proliferation of parameters, the ratio 
>^nnn/"^nn is kept fixed to 0.93, which is the value that 
results from ti = -1, = 1-3, ^3=^4 = -0.85^^ For 
the electron-doped case, the results were found to be not 
sensitive to changes in this ratio as long as the model re- 
mains in the regime with (tt, 0)/(0, tt) magnetic order 

The extended two-orbital Hubbard Hamiltonian is ex- 
actly investigated using the Lanczos algorithm2i2^ on a 
tilted -s/S X a/S cluster, as done in previous studies.— i^^i^'^ 
In spite of the small size of the cluster, this still requires 
substantial computational resources. More specifically, 
even exploiting the Hamiltonian symmetries the calcula- 
tion of the undoped-limit ground state of the eight-sites 
cluster still requires a basis with ~ 2-20 M states (slightly 
more demanding than a 16-site cluster one-band Hub- 
bard model), depending on the subspacc explored. Runs 
applying the Lanczos technique had to be performed for 
all the allowed momenta k of the cluster, and for all the 
quantum numbers under rotations and reflections (i.e. all 
the irreducible representations Aig, A2g, Big, B2g, and 
Eg of the symmetry group*^^), and also for all the z- 
axis total spin projections. In addition, the computation 
of binding energies for the case of hole doping requires 
calculations for a number of electrons N equal to 14, 15, 
and 16, varying U, Jh, and Jnn using a fine grid. For 
these reasons, the overall effort amounted to ^ 8, 000 di- 
agonalizations of the cluster, supplemented also by calcu- 
lations of dynamical properties, using a Penguin 128GB 
Altus 3600 computer. 



III. RESULTS 

It is important to remind the readers that the most 
commonly used models for the pnictides are usually de- 
fined in the so-called unfolded Brillouin zone .^i^'^'^^i^'^'^^ 
Due to the symmetry of the Fe-As planesji^i^ it is possi- 
ble to describe the pnictides using Fe-only effective mod- 
els where the As atoms merely provide a bridge for the 
electronic hopping between the irons. Under this approx- 
imation only one Fe atom is left per unit cell to describe 
these materials)^ As a result of these considerations, the 
number of orbitals to be considered is reduced by half, 
which is a computational advantage, and the size of the 
Brillouin zone (BZ) is doubled. For this reason, the mo- 
mentum in the unfolded zone is dubbed "pseudocrystal" 
momentumi^^ In order to relate the model results to ex- 
periments addressing the BZ corresponding to two Fe 
atoms, it is necessary to "fold" the extended BZ in such 



a way that the pseudocrystal momentum k = (tt, tt) is 
folded onto momentum (0,0). The physical difference 
between states with k ~ (0, 0) and (tt, tt) is that the 
first indicates a bonding and the second an antibonding 
combination of the d-orbitals in the two Fe atoms in the 
two-atoms unit cell. In the presentation of our results 
below, k will stand for pseudocrystal momentum. 



A. Phase diagram 

The relative symmetry between the undoped (A^ — 16) 
ground state (GS) and the A^ = 14 GS has been stud- 
ied with the Lanczos technique varying C//|ti| and Jn/U . 
The undoped GS was found to have momentum k = (0, 0) 
and it transforms according to the Aig representation of 
the Difi group, for all the values of Jh and U studied here, 
in agreement with previous results.— However, a surpris- 
ing result found in the present study of the hole-doped 
extended two-orbital model is the presence of many com- 
peting low-energy states not only with different symme- 
tries as in the electron doped case,— but also with differ- 
ent pseudocrystal momenta k. In other words, low lying 
states with both k = (0, 0) and (tt, tt) were found in our 
Lanczos investigation. This is compatible with previous 
mean-field approximation results that also reported low- 
energy spin-singlet pair states with momentum (7r,7r)i^ 

The competition among low-lying states with different 
symmetries and with different values of k is presented 
in Fig. [1] for the case Jn/U = 0.20, without the extra 
"J" terms. Numerically, it was found that the ground 
state for 14 electrons has crystal momentum (tt, tt). For 
small values of U this state is a triplet with A2g sym- 
metry (open circles in the figure). With increasing U, 
a transition (via a level crossing) occurs at U ~ 6|ti| 
to a spin-singlet ground state with B2g symmetry (open 
squares in the figure). However, it can be observed that 
there are states with k = (0, 0) that have very similar en- 
ergies. For example, for this pseudocrystal momentum, 
and in the weak coupling regime, a spin-singlet state with 
Big symmetry (represented with filled triangles in the fig- 
ure) is the closest in energy to the ground state, while for 
U > 3|fi| a spin-singlet state with Aig symmetry prevails 
(represented with filled diamonds in the figure). 

Similar results were found for all the values of U and 
Jh studied, i.e., the A^ = 14 ground state has total mo- 
mentum k = (tTjTt) but there are k = (0,0) states close 
in energy with a different symmetry. For this reason, the 
phase diagrams obtained by varying J-n/U and ?7/|ii| for 
both values of the pseudocrystal momentum will be pre- 
sented. 

The relative symmetry between the ground state with 
two electrons less than half filling with total pseudocrys- 
tal momentum k = (0, 0) and the undoped ground state 
is shown in Fig. [IK a), varying J-a/U and U/\ti\. A re- 
gion with symmetry -Big, indicated by the triangles, is 
found for small J7/|ti| (roughly J7/|ti| < 3) and moder- 
ate to large values of J-a/U . For larger values of J7/|ii|, 
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FIG. 1: (Color online) Difference between the energy of the 
lowest excited state with the symmetry and momentum in- 
dicated and the ground state. Full (open) symbols denote 
k = (0,0) (k = (tTjTt)). The results are obtained using the 
Lanczos algorithm for the two-orbital model in an eight-site 
cluster with 14 electrons (two holes doping), varying the Hub- 
bard repulsion U, and at a fixed Jh/U = 0.20. The inset 
shows a larger energy range in which the lowest lying state 
with each symmetry is displayed. The results shown in this 
figure are without the extra Jnn and Jnnn terms. 



the symmetry changes to Aig. A similar transition 
from Big to Aig (extended s-wave) has been found us- 
ing the RPA technique for an electron-doped five-orbital 
model at Jh = 0<^ The binding energy Eb, defined as 
Eb = E{U) + E{16) - 2E{15), where E{N) is the GS 
energy for N electrons, was also calculated. It was found 
that without the addition of Heisenberg terms there are 
no regions with binding. 



For the other case of a pscudocrystal momentum k = 
(7r,7r), the analogous numerical results are shown in 
Fig. [HI a). It was found that an A2g spin-triplet ground 
state, indicated by circles in the figure, dominates for 
large values of Ju/U and small C//|ti|. For the electron- 
doped model, an spin-triplet with momentum (0, 0) 
was similarly observed at large Jh and small U i^^^ For 
smaller Jh/U and larger C//|ti|, a spin-singlet ground 
state with B2g symmetry is the ground state. For this 
pseudocrystal momentum, the binding energy was cal- 
culated as well: binding was obtained for J^/U = 0.35 
where a spin-triplet ground state with symmetry A2g pre- 
vails (see open triangles in the figure). 
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FIG. 3: (Color online) Relative symmetry between the A^=16 
(undoped) and A'^=14 (with k = (vr,7r)) ground states vary- 
ing U and Jh/U. Circles denote spin-triplet states and 
squares i329-symmetric singlets, (a) Results for couplings 
Jnn ~ Jnnn = 0. Open triangles indicate binding, (b) 
Results for the lowest value of (Jnn, Jnnn) where binding 
appears with a fixed ratio Jnn/ Jnnn = 0.93.— 
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FIG. 2: (Color online) Relative symmetry between the 7V=16 
(undoped) and A'^=14 (with k — (0,0)) ground states, vary- 
ing U and Jh/U. Triangles denote Big-symmetric singlets, 
and diamonds Aig-symmetric singlets, (a) Results for cou- 
plings Jnn ~ Jnnn ~ 0. (b) Results for the lowest value 
of (Jnn, Jnnn) where binding appears with a fixed ratio 
Jnn/ Jnnn = 0.93.— 



B. Binding stabilization. 

To stabilize hole binding in the two-orbital model, 
we will proceed as in the previous investigation of the 
electron-doped caseii by adding extra Heisenberg terms, 
namely a NN coupling Jnn and a NNN coupling Jnnn as 
discussed in Sec. |IT1 As in the electron doped case, and 
as already explained, Jnn will be varied while Jnn / ^nnn 
will be kept fixed at the value 0.93iii 

The results for pseudocrystal momentum k = (0, 0) 
arc presented in Fig. [DJb), showing the symmetry of 
the hole-doped ground state for the lowest value of Jnn 
where binding of holes is achieved. The phase diagram 
remains largely unchanged by the addition of the Heisen- 
berg terms except for the Big region that has expanded 
slightly towards larger values of U. On the other hand, 
for states with momentum k = (tt, tt) the spin-triplet 
region virtually disappears (Fig. [3Jb)), except for those 
triplet states that already had Eb < at Jnn = 0, leav- 
ing behind a much larger B2g region in parameter space. 
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In Fig. SJa), the binding energy Eb vs. Jnn/U for 
states with momentum k = (0, 0) is shown for several 
values of U and at a fixed (realistic) Jn/U = 0.2. In- 
creasing Jnn eventually induces binding for all U^s. The 
value of Jnn/C^ where binding occurs decreases as U in- 
creases. Figure [5l[a) shows the same information but for 
the states with momentum k = (tt, tt), where a similar 
qualitative behavior is observed. 
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FIG. 5: (Color online) Results for states with total momen- 
tum k = (tt.tt). (a) i?B/|ti| vs. Jnn/U for different values 
of U/\ti\, at .Jh/U = 0.2 and Jnn/Jnnn = 0.93. (b) Phase 
diagram showing "Binding" and "No Binding" regions and 
the symmetry of the two-hole bound state varying (7/|fi| and 
Jnn/U, at a fixed Jh/U — 0.2. The shaded region indi- 
cates the "physical region" according to standard mean-field 
calculations.— The doted line is for Fig. [TJb). 



FIG. 4: (Color online) Results for states with total momen- 
tum k = (0,0). (a) EB/\ti\ vs. Jnn/U for different values 
of U/\ti\, at Jh/U = 0.2 and Jnn/Jnnn = 0.93. (b) Phase 
diagram showing "Binding" and "No Binding" regions and 
the symmetry of the two-hole bound state varying (7/|fi | and 
Jnn/U, at a fixed Jh/U — 0.2. The shaded area indicates 
the so-called "physical region" obtained from standard mean- 
field calculations that were compared with neutrons, trans- 
port, and photoemission experimental results.— The doted 
line is for Fig. [7] (a) . 



A study of the binding energy Eb and the relative 
symmetry between the N=16 and 14 GS's allows us to 
construct phase diagrams in the {U, J-^'n/U) plane. In 
Fig. mjb), typical results for the case J-r/U ~ 0.2 arc 
shownSl for the states with total momentum k = (0,0). 
The bound state has Aig symmetry in most of the bind- 
ing region, but a state with Big symmetry prevails at 
smaller U values (^ 3|<i|). In Fig. [5jb) the same infor- 
mation is displayed but for states with total momentum 
k = (tt, tt). In this case, the entire binding region, except 
for Jn/U > 0.3, has i?2g symmetry. All of the above 
symmetries appear inside the proper magnetic/metallic 
region of the undoped limit (indicated with shading in 
the figures) that were obtained in previous mean-field 
calculationsi^i^ extended to incorporate Jnnj^ 



C. Magnetism 

Since the two-orbital Hubbard model for the pnictides 
is not particle-hole symmetric, it is interesting to study 



how the nature of the doping, namely electrons vs. holes, 
affects the intensity of the magnetic order. In the actual 
materials, experimental results have shown that the in- 
plane resistivity of electron and hole-doped FeAs-based 
pnictides displays a larger anisotropy in the electron- 
doped easels Thus, it has been conjectured that the 
xz/yz magnetism is stronger in the electron-doped case, 
while in the hole-doped case it is weaker with a growing 
contribution of the xy orbital, disregarded in the two- 
orbital model, that forms the hole pocket around Mi^ 
A similar conclusion was reached via the FLEX approxi- 
mation for the case of electron and hole doping of a five- 
orbital Hubbard modeli^ 



The results for the two-orbital model studied here are 
shown in Fig. [6] where the magnetic structure factor S'(k) 
is shown in the undoped (crosses), electron doped (cir- 
cles), and hole-doped (diamonds) regimes, at fixed cou- 
plings J7 = 3, Jn/U — 0.2, and Jnn/C^ = 0.2, namely in 
the mean-field calculated "physical region" indicated in 
Fig. O While doping reduces the strength of the peak at 
k = (ti'jO), it is interesting to notice that the intensity 
is slightly more reduced in the hole-doped case. These 
results lend qualitative support to the notion that the 
magnetism in the xz and yz orbitals is stronger in the 
electron-doped case, and it becomes reduced when holes 
arc introduced. 
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FIG. 6: (Color online) Numerically calculated magnetic struc- 
ture factor 5(k), as a function of the momentum, using an 
eight sites cluster. Results for the undoped N=16, electron- 
doped N=18, and hole-doped N=14 cases are indicated, for 
couplings U/\ti \ = 3, Jh/U = 0.2, and Jnn/U = 0.2. 



D. Overlap Integrals 

In this subsection, the functional forms of the hole pair- 
ing operators that produce the hole bound states will be 
analyzed. With this goal, the overlap defined by 

(*(jv=i4)(k')|Ak'-k,«|*(iv=i6)(k = (0,0))) (3) 

was calculated using the Lanczos algorithm along the 
paths indicated by the dotted lines in the phase diagrams 
shown in panels (b) of Figs. 0] and \5\ Notice that for 
|\['i4(k')) the pseudocrystal momentum k' will take the 
values (0, 0) and (tt, tt) and, thus, a pairing operator with 
the appropriate k' — k has to be used to ensure a non- 
zero overlap. The ground state |\E'(iV)) in the subspace 
of N electrons was used, and the operator in Eq. (3) was 
defined as 



/(k)(a, 



lQ/3«k;,Q,t"k,/3,4,, 



(4) 



where (ik,a,(T destroys an electron with spin z-axis projec- 
tion cr, at orbital a = x^y, and with momentum k. The 
structure factor /(k) arises from the spatial location of 
the fermions forming the pair^^^ and ai arc the Pauli ma- 
trices {i = 1, 2, 3) or the 2x2 identity matrix ctq (* = 0). 
Note that ai and a2 imply an inter-orbital pairing. Over- 
laps for all the symmetries in Ref. [l^ and with NN and 
NNN locations for the electronic pairs, were numerically 
evaluated. 

In Fig. El^a), the overlaps for pairing operators with 
pseudocrystal momentum k = (0, 0) are presented for 
values of U and Jh along the dotted path in Fig. HJ^b). 
In the Aig region in Fig. Elja), we found that the same 
four pairing operators that have a finite overlap in the 
electron doped casoii also have one here. However, 



the relative strength of the overlaps differ. For consis- 
tency, we will use the same labeling for the operators as 
in Ref. [ill The Aig operator with the largest overlap 
is the operator (ii), i.e. the s± operator characterized 
by /(k)(Tj; = (cos fcj, cos fcj^)croi as in the electron doped 
cascr^ it is indicated by hollow diamonds in Fig. [TJ^a). 
However, in the hole-doped system the overlap for the 
pairing operator (iv) characterized by (cos — cos ky)a^ 
(hollow circles) follows in strength; this operator had the 
weakest overlap in the electron doped cascj^i The pair- 
ing operator (i) with (cosfca, + cosky)uQ (hollow squares) 
has an overlap almost as strong as in the electron-doped 
case. Finally, the overlap corresponding to the operator 
(iii)(sin fc^ sin fcy)(Ti (hollow triangles) is even more sup- 
pressed upon hole doping than upon electron doping. 

In the region where the pairs have Big symmetry there 
are three pairing operators with large overlaps : (viii) 
(cos fca; -|- cos fcj,)cr3 (solid circles); (ix) (cos fc^ cos fcj,)(T3 
(solid diamonds); and (x) (cos kx cosky)(7Q (solid 
squares). At small values of Jnn/?^, (ix) has the largest 
overlap amplitude followed by (x) and (viii). However, 
as J-mn/U increases (viii) overtakes (ix). 

For the case of pairing operators with pseudocrystal 
momentum k = (7r,7r), there is one contribution that 
clearly dominates, see Fig. [TJb): (vi') (cos fc^ cos fcj^)(Ti 
which corresponds to a NNN pair with i?2g symmetry. 
The prime in the label is used to remind the reader 
that the operator has a different pseudocrystal momen- 
tum from the B2g state with the same label discussed 
in the electron-doped case.— The only other nonzero 
pairing overlap occurs for (vii') (sin fcj; sin fcj,)(To and has 
a much smaller amplitude than (vi'). Interestingly, 
the nearest-neighbor B2g operator (v) characterized by 
(cos kx + cos ky)ai that had the strongest overlap in the 
electron doped caseii has zero overlap in the case studied 
in this manuscript. All the gaps for the pairing operators 
with B2c, symmetry have nodes along the x and y axes. 



E. Dynamical Pair Susceptibilities. 

To complete our analysis, the dynamical pair suscep- 
tibilities defined as 

/>oo 

P(c.)=/ dte-*(At (i)Ak,(0)), (5) 



were also studied in the state with A'^ = 16 for the pairing 
operators Ak.i introduced in Sec. IIII Dl Notice that the 
calculated spectral decomposition involves excited states 
with A^ = 14. The procedure described in Ref. [U in 
the context of the cuprates will be followed. As dis- 
cussed above, for A^ = 14 there are several low-energy 
states near the ground state that have different sym- 
metries. The dynamical pair susceptibilities show that 
most of these low-lying states have a large overlap with 
Ak,i|^'N=i6(0)) for Ak,i with the appropriate symme- 
try.' In Fig. M rcsuhs for U = 3\ti\, Jh/U = 0.2, and 
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FIG. 7: (Color online) Overlap {*(iV = 14) [ Ak,» |*(iV = 16)) 
vs. JNN/f/ for the indicated pairing operators, at f/ = 3 \ti\ 
and Jh/U — 0.2, for (a) states with total momentum k = 
(0,0) along the dotted path in Fig. |4jb), and (b) states with 
total momentum k = (tt, tt) along the dotted path in Fig.fS^b). 



Jnn/U = 0.10 are presented. Large overlaps with low- 
lying = 14 states are observed for operators (ii) and 
(iv) with Aig symmetry and (viii) and (x) with Big sym- 
metry, as well as for operator (vi') with B2g symmetry 
and pseudocrystal momentum k — (tTjTt). 

It is interesting to compare the results obtained for 
the dynamical pair susceptibility upon hole doping with 
those obtained for electron dopingiii In both cases, large 
susceptibilities for the low- lying states with Aig, Big, and 
B2g symmetries are found. This is remarkably different 
from the case of models for the cuprates where an analo- 
gous low-lying overlap analysis showed that dj.2_j^2-wave 
symmetry clearly dominates over all others.— 

Returning to pnictides, a similarity between the elec- 
tron and hole-doped cases is that the Big pairing op- 
erator that has the highest susceptibility, state (viii), 
is different from the Big pairing state for the cuprates. 
It corresponds to Cooper pairs mainly located on NNN 
sites, as opposed to the dominant NN contribution in 
the cuprates, and in the orbital basis used here the Big 
symmetry is realized by the orbital degree of freedom. 
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FIG. 8: (Color online) Dynamic pairing susceptibility for the 
pairing operators with total momentum k = (tt, tt) (operators 
with B2g symmetry) and with total momentum k — (0, 0) 
(operators with Big and Aig symmetry) (see text), at [/ = 
3.0 Jh/U = 0.2, and Jt^yi/U = 0.10. The vertical line 
indicates -Egs (14) - Egs{16)- 



In addition, the susceptibilities indicate that while NN 
pairs are favored in the electron doped case, NNN have 
larger susceptibilities upon hole doping for Aig and i?2g 
symmetries. 



IV. 



CONCLUSIONS 



The properties of a recently introduced two-orbital 
extended Hubbard model for the pnictides have been 
studied upon hole doping with the help of the Lanc- 
zos method. The results were contrasted with the previ- 
ously studied electron-doped case.— Due to the lack of 
particle-hole symmetry in the Hamiltonian, the results, 
as expected, are quantitatively different in both cases. 
However, an additional surprising characteristic of the 
hole-doped ground state is that it has pseudocrystal mo- 
mentum (7r,7r). In the reduced Brillouin zone represen- 
tation corresponding to the physical two Fc-atoms per 
unit cell description of the pnictides, having a nonzero 
pseudocrystal momentum means that the ground state is 
characterized by antibonding rather than bonding com- 
binations of the orbitals of the two Fe atoms in the 
unit cell. In terms of the pairing operators that are 
favored, it means that the pairs would arise from hole 
carriers located at the hole pockets at F and at M in 
the unfolded Brillouin zone. Interestingly, the five-orbital 
model for the pnictidesi^ shows that upon hole doping a 
hole pocket, absent in the electron-doped case, develops 
around M and the role of this pocket plays an important 
role in the properties of the hole-doped materials.— 
Our results may indicate that a simple toy model, such 
as the two-orbital model, could be used to study the role 
that a hole-pocket at M plays when multiorbital Hub- 
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bard models arc hole-doped. 

In spite of this difference in the pseudocrystal mo- 
mentum quantum number, there are several common- 
alities between the hole- and electron-doped two-orbital 
Hubbard models. The most important feature is that 
there are several low-lying states with different symme- 
tries close to the undoped ground state. For this reason, 
the symmetry of the doped states is strongly dependent 
on the actual values of the interaction parameters. Spin- 
singlet states that transform according to the irreducible 
representations Aig, Big, and i?2g were obtained both for 
hole and for electron doping. The richness of the phase 
diagrams unveiled here, and in the cited previous investi- 
gations, suggests that the symmetry of the pairing state 



in the pnictidcs is likely to depend on the material as well 
as on the type of doped carriers (electrons or holes) and 
on the density of dopants. 
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